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Abstract. We compute the field derivatives of a one- loop charge breaking effective potential 
and analyse their effect in its minimisation. The impact on charge breaking bounds on the 
CN . MSSM parameters is discussed. 

' Spontaneous gauge symmetry breaking occurs in the Minimal Supersymmetric Standard 

Model when the neutral - and obviously colourless - components of the Higgs doublets Hi and 
H2 acquire non zero vevs. The presence in the theory of many other scalar fields means there is 
^SJ ' no a priori reason why the minimum of the potential should be a charge and colour preserving 

■ one. If fields other than and H2 have vevs and for a particular combination of MSSM 

parameters the resulting minimum is deeper than the standard one we would therefore be in 
^ a situation where charge and/or colour symmetries were broken. This simple fact gives us, in 

Q_i' principle, a way of imposing bounds on the MSSM parameter space [Q]. There has been a great 

^ . deal of work done in this area most of it based on the analysis of the tree- level effective 

Q I potential along specific charge and/or colour breaking (CCB) directions. A fundamental point 

' in these works is that the CCB and MSSM potentials be compared at different renormalisation 

scales - this is based upon the work of Gamberini et al [p|], where it was showed the vevs 
^ , derived from the tree-level MSSM potential are a reasonable approximation to those obtained 

^ ' from the one-loop potential if one chooses a renormalisation scale of the order of the largest 

particle mass in the theory. Because that typical mass is different in the MSSM (of the order 
of tens or hundreds of GeV) and the CCB (of the TeV or tens of TeV order) cases, comparing 
both potentials at the tree-level order should in principle be done at two different scales. The 
authors of ref. [Q] determined CCB bounds including the one-loop contributions to the potential 
from the top-stop sectors, which they argued were the most significant ones. Recently the 
full one-loop potential for a particular CCB direction was calculated and used to restrict the 
parameter space of the MSSM. It was argued that comparing the MSSM and CCB potentials 
at different renormalisation scales neglected to take into account the field-independent part of 
those same potentials, vital to ensure their renormalisation group invariance |^]. The higgs 
and chargino contributions to the one-loop potential proved to be as important as the top-stop 
ones. This analysis was done at the typical CCB mass scale, so that, using the results of the 
computation of the one-loop derivatives of the CCB potential could be avoided. The results that 
were found had some renormalisation scale dependence, and it was then theorised that it would 
vanish if one performed the full one-loop minimisation of the CCB potential. In this letter we 
will undertake just that task. We rely heavily on the results of ref. H] and refer the reader to 
its conventions. We recall that we only consider the Yukawa couplings of the third generation, 
and the superpotential of the model is thus given by 



W = XtH^Qtn + XbHiQbR + KHiLtr + ^i/a^i 



(1) 



Supersymmetry is broken softly in the standard manner by the inclusion in the potential of 
explicit mass terms for the scalar partners and gauginos, and soft-breaking bilinear and trilinear 
terms proportional, via coefficients Ai and B, to their counterparts in the superpotential (^). 
At a renormalisation scale M, the one-loop contributions to the potential are given by 

where the Mq are the tree- level masses of each particle of spin Sa and Ua = (—1)^*°' (2so -|- 
1) Ca Qa- Ca is the number of colour degrees of freedom and Qa is 2 for charged particles, 1 
for chargeless ones. The "real" minimum occurs when the neutral components of Hi and H2 
acquire vevs vxjypi and V2/V^, the value of the tree- level potential then being 

^0 = \ {mlvl + mlvl) - BI2VIV2 + ^{9'^ + 9^) {vl - vff , (3) 

with mf = iTij^_^ + jj? and = m'jj^ + /i^. In the CCB direction we will consider the scalar 
fields tl and tji also have non-zero vevs, l/V^ and r/\/2 respectively, and the vacuum tree-level 
potential is now given by 

Vo = ^[vl{l^ + T^)+l^T^] - ^{ArVi + fiV2)lT + ^{mlvf + mlvl + mi + 

mlr') - B^,v,V2 + g (.i - - + 2 + {v'2 - vf + I'f . (4) 

The derivatives of this tree-level potential with respect to each of the vevs are very simple, but 
the same cannot be said for the one-loop derivatives, their total contribution given by 

Some of the squared masses' derivatives are trivial to calculate: that is the case of the top and 
bottom quarks, and the scalar partners of the second and first generation up and down quarks 
and electron and neutrinos (expressions (15) to (18) of ref. Q) - we present these results in 
the appendix. For the remaining sparticles the calculation is made more difficult by the masses 
being given by the eigenvalues of square matrices, sometimes as large as 6 x 6 - many of these 
matrices are much more complex than their MSSM counterparts due to the existence of charged 
vevs causing mixing of charged and neutral fields. For example, the "higgs scalars" of the CCB 
potential are in fact the result of the mixing between the neutral components of Hi, H2 and 
the fields tl and tr, their squared masses thus given by the eigenvalues of a 4 x 4 matrix. It is 
nevertheless possible to find analytical expressions for dM^/dvi, once the Ma themselves have 
been determined (which is easy to do numerically, where an analytical determination proves 
impossible) - this is accomplished by noticing that the particle masses are always given by the 
roots of an n*''-order polynomial (in our case, n = 2, 3, 4 and 6), 

F(A, Vi) = A" + ^A"-i + 5A"-2 + ... = , (6) 

with coefficients {A,B, . . . } generally depending on the vevs Vi = {vi,V2,1,t}. This equation 
implicitly defines the squared masses A in function of the Vi , so we have 



^ + ^A"-2 + 

dM^ _ dvi _ dvi dvi 



dvi dF_ nA"-i + (n-l)AA"-2 + (n-2)5A"-3 + ... 



■ (7) 

A=M2 



If n > 2 the final expressions depend on the numerical solving of eq. (P). For n = 2 it is possible 
to write down fully analytical expressions of the derivatives of the squared masses, but from a 
practical point of view it is better to use the recipe of eq. (|^). In the following we show how to 
calculate the coefficients {A, B,C, . . .} (and their derivatives) of eq. (^) for the several sparticles 
in terms of the elements {o, b,c, . . .} of their mass matrices. The derivatives of {a, b,c, . . .} are 
listed in the appendix. So, for a symmetric 2x2 mass matrix with diagonal elements a and c 
and off-diagonal element b, eq. ^ reduces to 

dc d{ac) 2b ^ 

^ dvj dvi dvi dvj ^^-^ 

dvi a + c — 2\ 

This is the case of the stop, sbottom and neutral gauge boson masses, the coefficients {a, 6, c} 
are given in eqs. (12), (14) and (21) of ref. Q, and their derivatives a simple calculation. The 
squared masses of the charginos are also determined by a quadratic equation [|, namely (from 
eq. (22) of ref. §), - A + B^± = 0, with 



^ 2 „. „, .. Ti/r\ , 1 „2 ,2 /^l 2 _2 , , 2\ , 1 2 „,2 , ^^2 \ 2 



^gl{l^V2 + M2Vl)lT . (9) 

All that remains is to calculate the derivatives of A-^± and B^± and substitute their values in 
eq. The charged Higgses are a mix between the charged components of Hi and H2 and the 
tau sneutrino, with a 3 x 3 mass matrix with elements {a±, b±, . . . , /±}, as shown in eq. (24) of 
ref. 1^]. The squared masses end up being determined by a cubic equation, 

- A±X^ - B±X - C± = , (10) 

with 

A± = a± + c± + f± 

B± = bl + dl + el - a± (c± + /±) - c± /± . (11) 

C± is, of course, the determinant of the mass matrix, but we end up not needing to calculate it, 
only its derivative. Adopting the convention X^, to indicate dX/dv, with v any of the vevs, we 
obtain 

A±,v = a±,-u + c±^v + f±^v 

B±,v = 2 {b± b±^y + d± d±^y + e± ± e±,„) - a± (c±,^ + f±^^) - c± (a±,„ + /±,„) - 
f± {a±,v + c±,^,) 

C±,v = 2{b±^^d±e± + d±^yb±e± + e±,„6±d±) + a±^yC±f± + c±^va±f± + 
f±,va±c± - d±{2d±^yC± + c±^vd±) - b±{2b±^yf± + f±,vb±) - 
e±(2e±,^,a± + a±,„e±) . (12) 

The squared masses of both the pseudoscalars and Higgs scalars are given by 4 x 4 matrices 
(eqs. (26) and (28) of ref. Q) with one set of coefficients {a, b, . . . j} for each case. The resulting 
fourth-order eigenvalue equation, 

X'^ - AX^ + B X'^ + C X + D = , (13) 



^One of the eigenvalues of the 5x5 chargino mass matrix is zero, corresponding to the r neutrino. This leaves 
a quartic equation in the masses, that reduces to a quadratic one in the squared masses. 



has coefficients 



A = a + c + h + j 

B = a{c + h + j) + c{h + j) + hj - - (f - - f - - 
C = a{g^ + f + - cj - ch - hj) + c{cf + + i"^ - hj) + 

h(b^ + e2 + g^) + j{b^ + + f) - 

2(/5? + bfd + dei + bge) . (14) 



Once again, for the calculation of the derivatives of eq. (^), we don't need the explicit value of 
D, the determinant of the mass matrix. With the same convention as before, we have 



B^v 



a{cv + + j^y) + c(a,t, + /^,^, + j,^,) + h{a^y + Cy + j^y) + 

j (a^y + Cy + h^y) " 2 (6 6,^, + dd^y + ee,t, + + 55,^, + ii^y) 

a^y {g^ + + i'^ — cj — ch — hj) + c^, (d^ + + — a j — a /i — /ij) + 

h^y {b^ + + g^ — a j — ac — cj) + j^y{b^ + d^ + p — ah — ac — ch) + 

2b^y {bh + bf - fd - ge) + 2d^ (dc + - bf - ei) + 

2e,„ (ec + e/i - di - bg) + 2f^y{fa + f j - gi - bd) + 

2g^y{ga + gh — fi — eb) + 2i^y{ia + ic — gf — ed) 



a,v (chj + 2fgi- jp 
h,y [acj + 2 beg — jb^ - 



- ci^ — hg'^) + c^y {ahj + 2 dei — jd^ 
ag^ - ce^) + j^y {ach + 2 bfd - hb"^ 



a% 



he') + 
-cd' 



+ 



2b ^y 


[bi" 


- bhj + fdj + geh 


- fei 


- gdi) + 


2d^y 


id 9' 


— cdj + f b j + cei 


- bgi 


- egf) + 


2e^y 


(ef 


— che + gbh + cdi 


- bfi 


- gdf) + 


2/,. 




- afj + dbj + agi 


— bei 


- gde) + 


'^9,v 


igd' 


— ahg + ebh + a f i 


— bdi 


- fde) + 


2 


iib' - 


- aci + a f g + dee - 


- bfe 


- bdg) 



(15) 



For the neutralinos, the sixth-order equation for the masses is made reasonably simple by the 
mass matrix (eq. (23) of |5i) having several zeroes. We thus have 



with 



-BvO A + CvO A + Z?vO A + -EvO X + -FvO A 

A, A, A. A, A, 







(16) 



B^o 



Ml + M2 

f (vl + /2 + r^) + i {g" + gl) (vf + vj + l') + g'\' + 



Ml M2 



-^{Mi +M2){vl + /2 + r2) 
g'^M2T' + ^{g'^ + gl) ^iVlV2 



1 



(/M2 + giMi){vi + vl + t') 



(3g'' - 52 - 2A^)t;i/T - //^ (Mi + M2) 



2^/2 



^{Mi + M2)vilr + A? { ^r2(.? + f+ r') + ^ (,'^ + gl) [vf + 

/2(^i - 2vf) + ^i(t;? + r^)] - iMiM2(.? + + r^) + ^ } + 

A. 
2^2 



(g^Mi - 3g'^M2)vi + (g'^ - gl)fiV2 



/2 



+ 



4 

A? 
V2 



2 

M1M2 



+ 9i)vl - M1M2 vilT + \l{ - ^(/M2 + [vf + + 



4 

(Ml + M2) vf - [M2 {vf +f + 



It 



/i Vl V2\ 



} 



+ 



Xr 
1 

- (5'' M2 + gi Ml)^l^l^ + ^ {gl viV2 - 2 ^i M2) 

-^(/M2 + g^2Mi)vilTvl + ^ [(/M2 + glMi)ii{f - vl)viV2 



2^?MiM2vl - ^V2T^ {4:flM2Vi + g2T^V2) 



+ 



A, 



*T ,2 



2V2' 



,,2 

/r, )\^ I 2 2 \ I /2 2 72 2 

{2flM2Vl + g2T V2)fJ-lT ^5 92^'^ 



(17) 



With such complex formulae, a check of the results is quite useful - because supersymmetry 
is softly broken, Str M'^ is a field-independent quantity, so we should have Str dM"^ /dvi = 0. 
In this manner we can check the mass matrices themselves and the consistency of our sign 
conventions. With the derivatives (^) computed, we can perform the one-loop minimisation 
of the CCB potential. We apply our formulae to the same simple model of ref. [^: one with 
universality of the soft parameters at the gauge unification scale and input parameters in the 
ranges 20 < Mg < 100 GeV, 10 < < 160 GeV, -600 < < 600 GeV and 2.5 < tan /3 < 
10.5, with both signs of the /i parameter considered. In the work of ref. [P, we had performed 
the tree-level minimisation of the CCB potential at a renormalisation scale M = g2\Ar\/ \t 
- which was shown to be of the order of the highest of the CCB masses - and, out of an initial 
parameter space of about 3200 "points" , CCB extrema had been found for roughly 40% of the 
cases. However, by repeating the process with the full one-loop derivatives of the CCB potential, 
we encounter drastically different results - only in almost 200 "points" does CCB seem to occur. 
These "points" are uniformly distributed according to the input parameters of Mq and mc, 
do not occur for values of //g close to zero (like in ref. ||5|) and occur mostly for tan/3 > 6 
and 150 < |74g| < 500 GeV. In figure (||) we see the reason for the discrepancy between these 
results and those of ref. [||: there we have plotted the evolution with the renormalisation scale 
M of the value of the one-loop MSSM potential ((Vq + AVi)^"'^"^'^^'^, with one-loop vevs) and 
the one-loop CCB potentials calculated with the tree-level CCB vevs {{Vq + AVi)^^^{vi) - for 
convenience, we divide it by a factor of 100) and the one-loop derived vevs {{Vo + AVi)'-''-''^ (vj)). 
To interpret this figure, we need to remember that Vq + AVi is not a one-loop renormalisation 
scale independent quantity Q, rather, in terms of the parameters Aj and fields (pj, the RGE 
invariant effective potential is given by 



y(M,Ai,0j) = niM,X,) + Vo{Xi,(t>j) + AVi{M,X,,cf>,) + o{h^ 



(18) 



^Verifying that Str is field-independent provides a check only on the diagonal elements of the mass matrices 
(except for the neutralinos) . But because we are using eq. (j^ to calculate the mass derivatives, this second check 
involves all the coefficients of the mass matrices. 



The only difference between the CCB and MSSM potentials is the different set of values for 
some of the fields which means the field-independent function Q is the same in both cases 0. 
Therefore, given that V is renormalisation scale invariant, we must have d(Vb + AVi)'^'^^/dM = 
d{Vo + /SVi)^^^^ /dM - this is certainly the case for the two one-loop minimised potentials of 
fig. (|^) (a plot of their renormalisation scale derivatives would show them to be almost identical 
for M > 1 TeV), as they run parallel to one another, but the one-loop potential calculated with 
the tree- level vevs is clearly diff^erent. It has a very strong M dependence, and the inequality 
(Vb + AVi)^^^{v^i) < {Vo + A^i)*^^^*^ is verified only for M < 4 TeV. Notice how, judging 
by the value of the one-loop minimised CCB potential, this "point" is not a CCB minimum. 
Unfortunately, we find that for those points that are identified as one-loop CCB minima the 
potential does not have the correct renormalisation scale dependence seen in fig. (|l|), as may 
be seen in fig. (^) - there, for a different choice of parameters, we obtain CCB potentials that, 
whether computed with tree- level or one-loop vevs, are strongly dependent on M. Although 
the one-loop vevs do seem to somewhat stabilize (Vq -|- AVi)*^'^^, one can expect its value will 
become greater than the MSSM potential for a higher renormalisation scale. Again, the finding 
of a CCB minimum becomes M-dependent. The reason for this seems to be the values of the vevs 
found - in the case of fig. (1) the one-loop vevs are smaller than 1 TeV, and very stable against 
variations in M. But, for the "points" where we find (Vq + AVi)'^'^^ {v}) < {Vq + AVi)^^^^^, 
the values of the vevs are much bigger than in the previous case and, as may be seen in fig. (^), 
change immensely with the renormalisation scale. In the same plot we also see the evolution 
of the tree-level vevs - remarkably they are rather stable with M, but the potential thereof 
resulting is still strongly M dependent. We must compare this figure to fig. (2.b) of the work of 
Gamberini et al whereas there, for a range of M of the order of the largest mass present in 
AVi, the tree-level and one-loop vevs coincide, in our CCB potential they simply touch in one 
particular point. We must add that for the seemingly perfect case of fig. (Q) the vevs do not 
even touch: Vi had a fairly stable value of about 4.6 TeV for the whole range of M, and also 
stable, was equal to ~ 0.8 TeV. 

In conclusion, we have shown that the one-loop contributions to the minimisation of a CCB 
potential have a large effect in both the values of the vevs and the potential itself. We found 
that the one-loop minimisation stabilizes both vevs and the potential against changes in the 
renormalisation scale, but only if the values of the vevs found are small ("small" in this case 
being inferior to about 1 TeV). In those cases, however, no CCB minima are found. We did 
find a small number of CCB minima, but the corresponding vevs had large values and the one- 
loop potential proved to be strongly M dependent. We believe the reason for this difference in 
behaviour is a breakdown in perturbation theory - for higher values of the fields, the one-loop 
contributions become too large and two-loop terms become necessary to achieve renormalisation 
scale invariance. As a result, the CCB minima found cannot be trusted - the corresponding 
value of the potential may well be smaller than the MSSM one, but perturbation theory is no 
longer valid and it is altogether possible the two-loop contributions would reverse that result. 
Also very important is the fact the tree-level derived vevs do not coincide with the one-loop 
ones for a range of renormalisation scale in the way described in ref. [^. The conclusion is 
that taking a renormalisation scale of the order of the largest mass present in AVi does not 
correctly reproduce the effect of the one-loop contributions to the potential, at least for this 
particular CCB potential and these values for the SUSY parameters - as the usual CCB analysis 
rely on this assumption, our findings cast doubt over their validity. Overall, the importance 
of performing a one-loop minimisation whilst studying CCB bounds cannot be underestimated, 
even if the results are not to our liking. To avoid the perturbative breakdown we encountered 
here we should study a CCB direction with lower typical masses, which suggests those directions 



^This was the argument used in ref. [q| to argue that both potentials should be compared at the same renor- 
malisation scale. 



associated with the top Yukawa couphng. We expect to approach this subject in the future. 



A Appendix 

We now Ust the non-zero derivatives with respect to the vevs {vi,V2, 1,t} of squared masses and 
elements of mass matrices. For the top and bottom quarks, the scalar quarks and leptons of the 
first and second generations and the charged gauge bosons, we have 



12 



Mi,,,, = - VI Ml^,^ = -o9 V2 Ml, = -g''l 



Ml,r = -\f r Ml^^^ = -1(5'^ + 3gj) v, M^^^ = ^(^'^ + 2>gl) v. 



M 



I , = -(-/ + 35i)^ Mj =-g'\ M? 



Ml,,, = --Agl - 5'') V, M? = -(gl - g'^) = -(gl + /) I 



Mir = -y^r M? = - M? = \g'^v. 



2' 



'ei,T 2^ 62, til 2^ 62,1)2 

Ml,, = -Y I Mir = 9" r Ml^^^ = \{g'' + gl) v, 

Ml,v, = -lig'' + gl) V2 Ml, = - gl) I Ml, = r 



M^,,, = ^glvi M^^^^ = ^glv2 M^, = ^gll (19) 



For the stop, sbottom and neutral gauge bosons, whose squared masses are given by symmetric 
2x2 matrices with diagonal elements a and c and off-diagonal element b, we have 

^(12 A| + / - 3gl) V2 ai, = -^(5" + 3gl) I 





= -12^' - 


3gl)vi 


«t>2 


%r 


1 r2 
= 7: Q T 
6^ 




Kv^ 


'^i,vi 


1 /2 






%r 


2 ,2 
= -3' " 






%,l 


12^^ 


3gl)l 


%r 


b,V2 


h 

= -7=2' 




bb,l 




= 1(6X1 -g")v. 


b,V2 


%r 


= 3' ^ 






boo, I 


= gl tan 


cos (2^^) / 





-{3\ - g )V2 ci, = -g I 

1(12 Xl - g'^ - 3gl) = ^(ff" + 3gl) V2 

— Q T Or = —= Ah 



- XhXrT bz^ = - 

2 '''^2 



1/2 1 ,2 

22 /2 

2^2 sin 9wl o,Q0.^ = 2g r 

5i 5i 



2cos2^w^''^ 2cos^ew''^ 



For the charged higgses, we have 





\i9" + d)vi 




-\{.9''^ - 91)^2 


a±,T = 


1(2X1 +9") r 


b±,Vl = 


1 2 

-^92'»2 




-\i9'^ -9i)vi 




1 , ,2 9\ 

-^K9 +92)V2 


C±,T = 


1 /2 

2^ " 


d±,vi = 




d±,r = 




^±,V2 — 




e±,T = 


-7=2' 


f±,Vl = 


\i9" + 9l)v. 


f±,l = 


li9" + 9l)l 


f±,T = 





For the pseudoscalars, 





= \i9"+9'2)v. 


"'H,V2 = 




+ 52) V2 




= 1(2X1 -g")r 


^H,vi = 


1, /2 


+ 92) vi 




= -l(9"-9l)l 


'^H,T = 


1 /2 

2 ^ 




^5,1 


" V2^^ 


fH,T = 


A. 






= \{4X'^+g" -gl)v. 


^H,V2 = 


1. /2 


- 92) V2 




= l(2Xl-g")r 


'^H,vi = 


V2 


T 






jH,V2 = 


1 ,2 

- g V2 




jH,r 











And for the higgs scalar s, 





l(9"+9l)v. 


aH,v2 


1 , .2 2\ 
= -^(9 +92)^2 


a_ff,T = 


1(2X1- 9") r 




= -\(9''^ + 9l)v2 


CH,vi = 


-l(9" + 9l)vi 


CH,V2 


= \(9'^ +92)V2 


CH,T = 


1 ,2 
2' " 


dH,vi 


= \(AXl+g''- 


dH,T = 


^A 
V2 ^ 


eH,vi 


= \(2Xl-g")r 


eH,T = 


1(2X1 -g")v. 


fH,V2 


= -\(9"-9l)l 


fH,T = 


Xr 

-7=2' 


9H,V2 


1 /2 
= 7:9 T 

2 ^ 



o 



(20) 

a±,l = \(9''^ +9^)1 
b±,v2 = ^92^1 

c±,i = ~(9'^ +gi)l 

d±,i = ~(2Xl - gl)vi 
9l 

e±,i = ^V2 
f±,V2 = -\(9''^ +92)V2 

(21) 

aH,i = \(^Xl + g'^ - gl)l 
Ch,V2 = \(9''^+92)v2 

V2 ^ 

Xr 

hH,i = \(9" + 9l)l 
■ _ ^- 

^H,V2-^^^ 

jH,i = \(2Xl-g")l 

(22) 

a„,i = \{AXl+g'^ -gl)l 

bH,v2 = -^(9'"^ + 92)n 

CH,i = -\(9'^ -92)1 

dH,i = \(^Xl+g'''-gl)v^ 

fH,i = -^(g''^ - gl)v2 

Xr 

9H,l = -^t^ 



9H,T = 


1 /2 


hH,vi = 




hH,V2 = 


-4(5 - 




hH,l = 


O 


hH,T = 


1(2X1- 9") r 


iH,vi = 






'^H,V2 — 


A. 
-7=2' 


iH,l = 


1(2X1- 9") r 


iH,T = 


i(2A? - 


9 


jH,vi = 


1(2X1 -g'')v. 


jH,V2 = 


2' 


jH,l = 


i(2A? - 


9") I 



>2 



jH,T = 3g T . 
For the charginos, from eq. (P), we have 

\'^,{vi,V2,l) = 92 (Vl , V2 , I) 



\±,T = X^T 



^X^,V2 



\9t(vl + l^)vi - ^glM2lT - gliiM2V2 



1 



A? 



1^92 vlv2 + y32t'^^2 - ^giulT - gl^M2Vl 



B^±,i = gl (;u2 + ^vi^ I - ^9i(f^V2 + M2V^)t 



r2 , 92 „,2 



A. 



K Mi + ^7^2 - -^gi(^V2 + M2Vi)l . 



V2- 



Finally, for the neutralinos, from eq. ([T7|), we find 



Cx°,V2 



^X°,l 



'^X°,r 



^X°,vi 



1(2X1 + g'' + gj)v. 
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Figure 1: Evolution of {Vq + AVif^^{vj), {Vq + AT4)^^-^(v°)/100 and {Vq + AVi)^^^^ 
(solid, dashed and dotted lines respectively) with the renormalisation scale M. Notice how the 
two one-loop minimised potentials run almost parallel to one another. For this case, Mq = 80 
GeV, mo = 40 GeV, Ag = 400 GeV, tan/3 = 6.5 and no = -1-436. 
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Figure 2: Evolution of (Fq + AVi)^'^^{vl), {Vq + AFi)^^^(uO) and (^o + A^i)^^^^^^ (solid, 
dashed and dotted lines respectively) with the renormalisation scale M. The strong variation 
with M of the vevs displayed in fig. causes the similar dependence of the CCB potential 
here. Mq = 100 GeV, mo = 160 GeV, Ac = -600 GeV, tan/? = 10.5 and hg = -2.9823. 




Figure 3: Evolution of the vevs vi (from one-loop - solid line - and tree-level - dashed line - 
minimisation) and V2 (from one-loop - dot-dashed line - and tree-level - dotted line - minimisa- 
tion) with the renormalisation scale M, for the same choice of parameters of fig. (^). Notice the 
large difference in value between the tree-level and one-loop vevs. For very high values of the 
renormalisation scale, the one-loop vevs begin to stabilise. 



